THERMAL GAUGE THEORY, EFFECTIVE ACTION AND Z(N) 
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Abstract 



One of the salient features of high temperature gluodynamics is the Z(N) groundstate 
degeneracy, key to the understanding of the transition. On the other hand impressive 
progress has been obtained with effective theories, that seemingly lack Z(N) symmetry. 
We give a simple unified description of the two. As a very useful byproduct we get a 
qq , natural definition for the Debye screening mass. We show its expediency by calculating 

' the next to leading order in the Debye mass. The sign of this correction gets a simple 

**"~^ explanation. 



1 Introduction. 

Thermal gauge theory is the main tool to study equilibrium properties of plasmas in 
elementary particle physics. Their relevance for ongoing experiments at CERN, for the 
physics at RHIC and LHC motivates their study, both through analytic and numerical 
methods. 

Iruracetit years very useful results have been obtained through effective thermal the- 
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■ ories u'lji I. Properties of the bulk plasma phase have been elucidated, leading to new 

insights. Markedly absent in effective theories are the Z(N) symmetries, that gave so much 
impetus to the understanding of the deconfinisg phase transition □. An inroad is made to 
the relation between these two aspects in ref. 13. 

The present note gives a simple unified picture, and some concrete consequences of this 
picture in perturbation theory fl 

We concentrate on pure SU(N) gluodynamics at some temperature T, well above the 
critical temperature T c , where deconfinement takes place. Such a theory is characterized 
by the value of A„i„ and T. Their knowledge determines the coupling g. In perturbation 
theory at high T (so g = g(T) small) three different scales come into play: T, gT , and 
g 2 T. The electric screening length mr> is of order gT, and there is a magnetic glueball 
mass of order g 2 T that limits the applicability of perturbation theory. 

Now we consider a 3D box of size L 2 r L z , elongated in the z-direction, containing the 
plasma. Its statistical properties are described by the Gibbs sum over physical states: 

cxp -F/T = V cxp -H/T (1) 

phys 



a It is known for some time that convergence is slow for certain observablesQOQ. So its use for phenomenology 
is not clear. Nevertheless we feel that control of these aspects still warrants efforts. 
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Physical means we impose Gauss' law everywhere inside the box. Only at the sides of 
our box we may drop Gauss' law and impose special boundary conditions, specified below. 
The Hamiltonian equals 

H = J dxTr (g 2 E 2 + -^B 2 ^j (2) 

The sum can be rewritten as a Euclidean path integral with period 1/T in the time 
direction: 

exp-F/T = J DAoDAexp--^S(A) (3) 
The Z(N) invariance is naturally defined in terms of the Wilson line: 

f f 1/T 
P(Ao) = -TrPexpi J A (r,x) (4) 

Now the allowed gauge transformations are those that leave the action invariant. Apart 
from periodic ones, we have transformations that are periodic modulo the centergroup 
elements exp ik^-. Those are the only ones leaving the action invariant. So the phases 
C"(af) of the eigenvalues of the loop are gauge invariant modulo multiples of 2n/N for all 
i running from 1 to N. What is important is that any adjoint field- or more generally 
any representation with zero N-ality- stays periodic (or anti-periodic) under any allowed 
transformatior|| 

Anticipating on the wall profile p(z) in the z-direction we define the phase average 
C l {z) = y±zJ (fx±C i (x). 

The free energy of our system ( ^ can be written as: 

exp-F/T = J Dp{z) J DA DA S{p(z) - C) exp ~^S(A) (5) 

So we split the path integration into an integration over profiles p l and the remaining 
variables (for readability's sake we left out the color indices in eq. ^. The remainder is an 
integration over all potentials with the gauge invariant delta function constraint. For very 
large size of the box we can write it as: 



exp-|L? r L z - 1 S ' f ^ L%. = [ DAoDA S(p(z) - C)exp~S[A) (6) 
T g 2 N T J g 2 

The first term in the exponent is the free energy per unit volume and does not depend 
on the profile p, which describes surface effects. 

The second term is the center of our interest. It contains the information on the 
dynamics of the wall. In the limit of very large box size Lt r becomes so large that only the 
profile with minimal action will contribute. So in that limit knowledge of the minimum 
of S e ff is sufficient. Note that the effective action takes the same value for all p that 
differ by a centergroup transformation, because a gauge transformation periodic modulo 
a centergroup element will not change the value of the integral in eq. [| Only the phase of 
the Wilson line in the argument of the delta function will change by a shift k^S-. Boundary 
conditions on p at the edge of the box sk, z = ±kf- fix the minimum or domainwall as in 
fig. 1 and correspond to fixing the value of the Wilson line to be 1, respectively exp i 2 ^-. 

A last remark concerns the physical meaning of the profile, where it merges into the 
groundstate of the plasma, on the far left or rig ht in fig. |l|. From its definition in eq. a 

b Also for non-zero triality representation the action stays invariant. But the periodicity is affected 
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it is a gauge invariant version of the Aq potential. In the Abelian case it is just the line 
integral of Aq, and we know that it decays exponentially in the plasma. Hence we will take 



(7) 

as the definition of the Debye mass in the non-Abelian case. This definition is gauge 
invariant and non-pjerturbative. It is odd under Euclidean time reversal, so cannot excite 
magnetic glueballs u. We will see that it coincides in perturbation theory with the usual 
value of the Debye mass as defined.^ the two point functions or by the correlation of the 
imaginary part of the Wilson lineut-j. 

2 Computing the effective action 

In this section the strategy of our computation is laid out. We suppose the coupling g 
is very small. Apart from the three scales in the theory: T, gT and g 2 T there is a scale 
due to the profile p, which is sliding. As we saw in the previous section, Z(N) invariance 
introduces a periodicity in tte profile, so p/T = O(l) or less. 

We follow the idea in ref.tj and integrate out the heavy degrees of freedom first. Heavy 
means 0(T). Then we are left with an effective action which is at most three dimensional, 
because the heavy degrees of freedom include all non-static modes. 

We can distinguish two cases, depending on the value of p/T. 

(i) p/T — O(l). Then we can integrate out all degrees of freedom and find S e ff as a 
function of the profile p directly. 

(ii) p/T = O(g) or smaller. All non-static degrees of freedom are integrated out, but 
not the static ones. The resulting 3D action containing the static variables has a 
built in scale, the Debye mass. It is an action with non-local vertices. Only when we 
are looking at distances much larger than the Debye scale we can approximate the 
vertices by local ones. Hence we will treat the case p/T = O(g) later, and discuss 
only the 0(g 2 ) case in this letter. 

3 Perturbative saddle point and integration of heavy modes 

In order to integrate out the heavy degrees of freedom we need perturbative information 
on the effective action, so we have to search for the saddle point. This saddle point cannot 
be anything else than the a priori profile p appearing in the constraint. A simple analysis 
of the constraint shows this immediately, see below eq. |^. This means the fluctuations Q 
around the profile can be written as A^ = p<5 M ,o + gQn, where p is a diagonal traceless 
N x N matrix, and Q any Hermitian traceless matrix. Then we have to fix the gauge, and 
the obvious choice is background gauge : 

jTr(D(p)^Q^) 2 . (8) 

To deal effectively with the z-dependence of our background field p(z) we will do a gradient 
expansion. This is justified, because the classical action of our profile is proportional to 
1/ g 2 (dzp) 2 ■ The next term is of quantum origin (the log of the fluctuation determinant) 
and is therefore O(l). Hence the minimum profile will have a gradient of order g to balance 
the two terms. 

So we will split the profile p(z) — p + 8p. The first term is constant. The second term 
is expanded out and gives rise to vertices. Were it not for the constraint, tha resulting 
Feynman rules would be the usual background field gauge rules (see e.g. ref. u) with all 
time derivatives do replaced by their covariant counterparts Do(p), since p can be of the 
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same order as do, that is, 0(T). So propagators for off-diagonal quantum fields change 
accordingly, and specially for the constant modes the profile acts like an infrared cut-off. 
Diagonal quantum fields do not feel the difference between the two derivatives, since p is 
diagonal, so commutes. So for SU(N) only N-l massless excitations remain. 

To see the consequence of the constraint in more detail we first expand the constraint 
in terms of the coupling g. 

After some algebra we get 



)-/'''•: —J <!■'■ [ ~Qa(*) + |L \^^(Ql J &n)Ql\$, -n)j + <>:<,■:■ j 

(9) 

in the traditional notation Q a = TrQX a , Q lJ = TrQ\ l \ with \ a the Gell-Mann matrices 
and A lJ the matrix with only one non-zero entry in the (ij)th place, all normalized to 1/2. 
The presence of the momentum component gfj(p) = 2n Tn+p 1 —p? reflects the non-locality 
in time of the constraint. 

So to lowest order the constraint tells us not to integrate over fluctuations that are 
diagonal and constant in the transverse directions. This is to be expected, because they 
are zero modes of the profile p(z). The only linear terms in the classical action in eq ^ are 
these same zero modes, so we have justified the profile p(z) as the saddle point. 

Now we expand the 0(g 2 ) terms out of the delta function. To reinstore the delta 
function constraint we have to do a partial integration with respect to the zero modes. 
This is going to generate couplings of the 0(g 2 ) term of the constraint to all derivatives of 
the action. These couplings stem from the non-linearity, i.e. the gauge invariance of the 
constraint and are therefore instrumental for the consistency of the approach. 

In fig. 2 we show the lowest order couplings. The b couplings are renormalizing ,±he 
kinetic term and the q terms the potential terms in S e fj much in the vain of Belyaevu. 

4 Results of integrating heavy modes 

We now integrate the right hand side of eq. ^ over the non-static modes Q. The result will 
be of the form: 

/ i s eff ( P ) r2 \ f nnfn , f s (3) (Q(Q),p) i s eff 2 \ 

ex P 7T7 L tr = / -DQ(0)exp —LLL tr (10) 



g 2 N T " J J \ T g 2 N T 

In the limit of small external momenta the three dimensional action can be written 
in a local form: 

5" (3) (Q,p)= j dx l/2j2 Tr (D(Q)kQi) 2 + Tr(D(Q k )Q ) 2 + Tr ([ip, Q M ] 2 ) 

+ jTr(d k Q k ) 2 +m 2 D {p)TrQ 2 + 8S (3) (11) 

With the profile p set to zero it reducesi-Jjenthe effective action used in computations 
concerning observables in the groundstate j].ali3 

The coefficients in S (3) are of order g 2 (Debye mass) or g 4 (SS^ 3 ') and stay so when 
p is turned onlij. For the next to leading order in the Debye mass we will only need the 
coefficient m 2 D (0) — ^g 2 T 2 . The remaining terms in eq. [id are the truncated kinetic and 
potential terms in S e ff(p). 



c Usually the gauge fixing term in the static sector s«k4nhence ghost terms are dropped: the static sector is 
considered as a non-perturbative sector. Here, like inoHj we keep the gauge fixing to pull out precisely the 
interesting-and still available-perturbative information from the static sector. 
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4-1 Profile of order T 

The p-induced mass terms are eventually getting larger than the Debye mass: in case the 
profile p/T is of order one we can integrate out the static modes on the same basis as the 
non-static modes, and change S e ff(p) to S e ff = K + V. K is the kinetic term and V the 
potential. Since these has been analyzed including twO|4cjop-|Order we shall be brief and 
refer the reader for more details to the original papers OLftJ. One obtains for all those 
values of p a perfectly Z(N) invariant profile (see fig. |l|). 




Figure 1: Profile of the wall. Dashed horizontal lines separate the q = O(l) region from q = O(g). 

We give the explicit expression for SU(N). We limit ourselves to the valley of minimal 
action. It is parameterized by one parameter q, that determines the matrix p as: 

p = 2nTq diag(l/N,l/N,...., (1 -N)/N) (12) 

The potential part V reads: 

V(q) = (N- l)(2nT) 2 m D (0) 2 (l - 5^^) J dzq\\ - q) 2 (13) 

with mo(0) the perturbative one loop Debye mass in eq. [y] and q mod 1. The first term 
comes from the fluctuation determinant. 

The second term comes from the two loop free energy graphs and from the renormal- 
ization of the Wilsonline through the q- vertex in fig. [| The lattep, is a typical example of 
how the vertices from the constraint restore gauge independenclffl. 



The result for the kinetic term reads: 



K= (N-1)(2tyT) 2 dz[l + 



11 g 2 (T)N 

3 (4tt) 2 



(</>(<?) +V(l-?) + l) (d z q) 



(14) 



The first term is the classical term. Combining this term with the lowest order term in 
V we find indeed that for small values of the profile the perturbative Debye mass controls 
its behavior for large z. So our definition of the Debye mass at the end of section [l] is 
justified. 
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Figure 2: Lowest order additional vertices coming from the constraint. The lines attached to the empty square 
are the Q's in eq. ^. The lines attached to the black square are due to differentiation of the three- vertex of the 
invariant action as explained below eq. N. Vertex q is 0(g 2 ), vertex b O(l). 



The second term is obtained from the graphs in fig. [j As is well known, these graphs 
give at zero temperature a p-independent result (at zero temperature we deal with an 
integration instead of summation over frequencies, so all p dependence is absorbed). All 
what stays is coupling constant renormalization. 

At finite temperature we are left with q dependence through the digamma functions. 
Gauge dependence in these finite parts is absorbed by the b-vertex induced contribution 
from fig. g. Since we integrated out the fluctuations at scale T the typical coupling constant 
renormalization coefficients are expected. The renormalized coupling in the results is 
obtained by dimensional regularization, and denned by: 

Some comments are in order. The results for K and V are the leading terms in a 
gradient expansion, where the gradient is units of the temperature T or the profile p. 

Trying to obtain the next to leading order for the Debye mass from eq. leads to 
disaster, since the digamma function behaves like ^ at small q. It is due to the static 
contribution, and an artifact of our use of Feynman rules not including the effective 3D 
action with the Debye mass. The latter is only justified when q = 0(1). 

4-2 The profile forp = 0(g 2 ), and corrections to the Debye mass 

As we just saw the kinetic and potential terms are in this region dominated by the static 
term in K. We know from Linde's vetocj that we cannot take the cut-off p/T smaller than 
g 2 , otherwise perturbation theory makes no sense anymore. From the explicit form for K 
we see that we will get the leading contribution for these values of p. So we computed 
the first four graphs given in fig. with the rules from S' 3 ' . In this action propagators 
involving Qo have now the Debye mass dominating the mass term induced by the profile. 
In the propagators where the Debye mass is absent the profile still provides a cut-off 
2nTq = m s . This renders our approach unambiguous. We have calculated with these 
rules the four diagrams in fig. ^ and for any value of the gauge parameter £. The leading 
result comes entirely from the first graph in fig. ^| It is for all £ equal to: 

277V 2 — lo s 

(47T) ir 



+ m s + mD(0) 



-ir + m s + mo(0) 



(16) 



The logarithmic term gives the cuts starting at ±i(m s + m£>(0)) in the momentum r 
of the incoming line. Since the momentum r equals imo(O) the contribution to the kinetic 
term becomes: 

K=(N- l)(2nTf f 1 - 2^^- log f ^)) (8 zq f (17) 
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Figure 3: 1-loop corrections to the kinetic term. Dashed lines are ghost lines, straight lines Q-lines, wavy lines 
background fields. The last graph on the right comes from the vertex b in fig. El 

Plugging this result into the equation of motion we get as correction to the Debye mass 

rno=rn D ( Q )(l + ^^loJ^m) (18) 
\_ 47t mn{0) \ rn s J J 

To leading order the argument of the log can taken to be 1/g. This is the result of Rebhan§. 
5 Conclusions 

We have given a unified treatment of the Z(N) invariant action and the effective action. In 
effective action parlance: in between two Z(N) vacua constant modes get heavy too, and 
can be integrated out with out recourse to non-perturbative physics. Near the Z(N) vacua 
the constant modes get light, and if they are light enough- 0(g 2 T) (but not lighter)- one 
can use the 3D effective action to compute the perturbative effects. 

Our approach brought two boons: a manifestly gauge invariant and non-perturbative 
definition of the Debye mass, and a gauge invariant infrared regulator. It warrants the 
computation of higher order perturbative effects in the surface tension, and in the bulk 
free energy at 0(g 6 ). 

The sign of the correction in eq. ^ is interesting. It is correlated to the sign of the 
coupling constant renormalization. At scales T it is opposite in Abelian and non-Abelian 
theories. This sign difference still persists at scales g 2 T, where we evaluate the Debye 
mass. Only the magnitude changes, both numerically and in the order of the coupling. 
This correlation is manifest in our choice of observable. 

An important issue is whether the <?log(l/(?) correction is not reproduced in higher 
order, admitting higher order vertices in the effective action S 1 ' 3 '. That is, one hopes to 
have at most O(g) corrections. In the context of the two-point function definition one has 
never found any counter examples. With the definition advocated above it may be easier 
to find a truly satisfying argument. While the non oerturbative O(g) corrections dominate 
the g log(l/g) corrections till very high temperatureu, the center of the domainwall remains 
accessible to perturbation theory at relatively low T. 

It is an open question whether there exists a version of our definition of the Debye mass 
that would permit us to study on a 3D lattice the wings of the domainwall. Breaking Z(N) 
symmetry by adding a complex representation does not affect our method of calculating 
bulk properties. 
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